term=4; n=400;

data(0.1);

A=matrixA(term,n);

b=source(term,n,0,1);

u=A\b;

clear A b

maxL2err(u,n,term);